Positive Magnetoresistance of Composite Fermions in Laterally Modulated Structures 
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Adopting the mean-field composite fermion picture, we describe the magneto-transport properties 
of a two-dimensional electron gas with laterally modulated density around filling factor v = 1/2. The 
occurrence of a strong positive magnetoresistance at low effective magnetic fields as well as Weiss 
oscillations, which were observed in recent experiments in such systems, can be explained within a 
semi-classical Boltzmann equation approach, provided one goes beyond a second order approximation 
in the modulation strength. 

PACS numbers: 71.10.Pm,73.40.Hm,73.20.Dx 

Transport properties of strongly correlated two- 
dimensional electron systems (2DES) in high magnetic 
fields, where the lowest Landau level is about half filled, 
seem to be surprisingly well described by the "Compos- 
ite Fermion" (CF) picture. This model is motivated by 
a singular Chern-Simons gauge-field transformation that 
maps the electron system near filling factor v = 1/2 onto 
a metal of quasi-particles (the CFs) . They can be inter- 
preted as electrons to which two magnetic flux quanta 
have been attached Q. In the usual mean- field approxi- 
mation, the CFs are considered as having the same den- 
sity ncF — n e \ = n as the electrons and as moving 
in the weak effective magnetic field B e g = B — 2$o™> 
where B is the external applied magnetic flux density 
and <f>o = h/e the flux quantum. Halperin, Lee, and 
Read H predicted that, at low temperatures, these CFs 
fill in momentum space the spin-polarized states within 
a well defined Fermi circle of radius fcp = v47m. There- 
fore, one expects that transport properties, which are 
governed by elastic scattering at the Fermi energy, can 
be calculated without considering explicitly the mutual 
interaction between the CFs, by analogy with the Landau 
theory of the electron Fermi liquid in zero and weak mag- 
netic fields. Indeed, interesting features of the magneto- 
resistance of 2DESs (near v — 1/2) in lateral superlat- 
tices, such as two-dimensional anti-dot lattices or one- 
dimensional (ID) superlattices, created either dynami- 
cally by the application of surface acoustic waves or stat- 
ically by surface etching, have been reproduced astonish- 
ingly well by calculations describing the CFs as classi- 
cal, non-interacting particles moving in suitable effective 
fields @-||. 

Here, we focus on static ID density modulations 
and compare experimental dc-transport measurements 
M with calculations based on the linearized Boltzmann 
equation (LBE) . The purpose of the present letter is to 
demonstrate that all characteristic features of the mag- 
netoresistance curve near B — B x / 2 = 2<£>o™ can be re- 
produced in the quasi-classical CF picture, provided one 
incorporates all modulating fields and solves the LBE be- 
yond the Beenakker-type approximation (B A) M , which 
allows an analytical solution and has been employed in 



previous work ^,g|Jg]. These features are: (i) a pro- 
nounced V-shaped minimum near B 1 i 2 , explained by 
"channeled orbits" which are omitted in the BA, (ii) 
shoulders or minima related to commensurability effects 
(Weiss oscillations due to drifting cyclotron orbits), (iii) 
a steep increase of the resistance with |-B e ff| between the 
structures caused by commensurability effects and the 
oscillatory structures due to the fractional quantum Hall 
effect, and (iv) an asymmetry in the slope and magnitude 
of these steep resistance flanks with respect to B e g = 0, 
due to interference effects between the direct electrostatic 
density modulation and the induced modulation of -B c ff ■ 
The material parameters needed to achieve agreement 
with the experiment around v = 1/2 are consistent with 
those needed to explain the commensurability effects at 
low magnetic fields [B < 0.5 T). 

Recent work by Mirlin et al. M , based on an analytical 
solution of the LBE within the BA for a special model of 
anisotropic scattering B , obtained reasonable results for 
the features (ii) and (iii), but not for (i) and (iv). Their 
approximation misses relevant physics at low B c g, where 
the CF picture is expected to apply best, and gives a 
rather poor fit to the experiment for very small |B e ff| 
(see Fig [j] below). 

To describe the resistance for B < 0.5 T, we follow 
Ref. U for pure electric modulation, and treat the 2DES 
as a degenerate Fermi gas of non-interacting particles, 
with charge — e and effective mass m*j, average den- 
sity n, Fermi energy Ef — n/Df, and density of states 
Dq 1 = m*i/(7r7l ). The etching of grooves into the surface 
of GaAs-heterostructures is assumed to produce an exter- 
nal electrostatic potential energy T^ oxt (x) = Vq** cos qx of 
period a — 2ir/q in the plane of the 2DES. It is screened 
by the 2DES and this leads within the Thomas-Fermi 
approximation to a potential energy V cl (x) = V^cosga: 
with Vo cl = V oxt /[l + 2/(<jaJj)], where a* B = l/[Ke 2 Df/K] 
is the effective Bohr radius. For GaAs, k = 12.4 and 
a* B rs 10 nm <C a ~ 400 nm. Consequently, the relative 
modulation strength is e c i = V^/Ef ps V r oxt Ka/(e 2 n), 
and the modulated electron density is n(x) = n[l — 
e e \cosqx], independent of m*j. 

To describe the resistance of CFs, we modify the ap- 



proach of Ref. M. We treat the CF system as a de- 
generate Fermi gas of non-interacting particles, with 



charge — e, effective mass m 
Fermi energy E F 



CF 
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CF _ 



cf O' average density n, 
h/DQ F , and density of states 
iQ ¥ /{2-n:fi l ), which obeys Newton's equation, 
ra^ F v = — e[F c ff + v x (B e se z )], with effective electric 
and magnetic fields. We use the LBE to calculate the re- 
sponse to an external homogeneous electric field E^ ). In 
the absence of E^ ), F c ff = W(x)/e is determined by the 
screened modulation potential, which we parametrize as 
V(x) = ecFEp F cos(qx). Since the static Thomas-Fermi 
screening of the CFs should be equivalent to that of elec- 
trons (with Dq F > Dq 1 instead of Df), we expect ccf ~ 
e e i, so that the modulated density n(x) — n[l— ecF cosqx] 
for B near -B1/2 is the same as that for small B within 
this approximation. This yields the effective magnetic 
field B Q g(x) = B — 2& n(x) = B c e + SB c e(x), where 
6B c g(x) — 6cf -B1/2 cosqx is the magnetic modulation. 
The linear response to E^ ) is carried by CFs at the 



Fermi edge, (m* CF /2)-v 2 + V(x) = Eg* . Due to the trans- 



= E§ F . 

lational invariance in y direction, the (suitably scaled) 
zero-temperature distribution function $>(x, p) 0] de- 
pends only on two variables, the position x and the po- 
lar angle p of the velocity v(x, <p) — v (x) (cos ip, simp), 
with v(x) = v F [l - V(x)/E§ F ] 1 / 2 , v F = hk$ F /m* CF , and 
fcp F = [Airn] 1 / 2 . We normalize $(x, <p) so that the cur- 
rent density j(x) reduces to 



m 



e 2 DZ F 



■^■v(x,tp)$(x,<p). (1) 



A current of CFs implies motion of flux tubes and pro- 
duces the Chern-Simon electric field 10 



E ob = -(2Ve 2 )jxe 2 , 
which, in the linearized Boltzmann equation, 

r>$ - C[$] = v(x, p) (e (0) + E cs (x) 



(2) 



(3) 



adds to the external driving field E^ ) . Therefore a self- 
consistent solution of Eqs. ([l])-([|) is necessary. The drift 
operator 

V = v(x)coscpd x + (u> c + (j m (x) +u> e (x,<p))d ip , (4) 

contains the average value and the periodic component 
of the effective magnetic field via lu c — eB c ^/m* CFl and 
u> m (x) = eSB c g(x)/rriQ F , respectively. The electric mod- 
ulation enters in ui c (x,p) — —smtpdv/dx and v(x) M. 
The collision operator is taken to be the same as for the 
unmodulated system and is written in the form 



cm 






P(p>'-p)l<P(x,p')-<I>(x,p)}, (5) 



where the differential cross section is parametrized by a 
relaxation time r and a dimensionless kernel P(p)- For 



actual calculations we take P{<p) 
< b < 1 and 

P p (p) = [(2Pp\) 2 /(2p)\]cos 



b+(l-b)P p (ip) with 



2l> 



fa/2) 



(6) 



which has the finite Fourier expansion P{<p) = 
J2n=o7n cos rnp with 70 = 1, and j n = 2(l-b)(p\) 2 /((p+ 
n)\(p—n)\) for n > 1. For b — 1 andp = 0, P(p) describes 
isotropic scattering. For b < 1 the fraction (1 — b) of the 
total scattering cross section is due to anisotropic scat- 
tering. With increasing p, P p (p) is increasingly stronger 
peaked in forward direction. 

To solve Eqs. (|l|)-(@), we calculate the distribution 
function of the homogeneous, unmodulated CF system 
first. This yields the resistivity tensor g h with compo- 
nents given by g\ x = g\ y = q$ f and g xy = -g\ x = 
wtotTtrgo F ; where g$ F = m* CF / (e 2 nT tr ) and r tr = r/(l - 
7i/2) is the relevant CF transport scattering time. Al- 
though the CFs move in the effective magnetic field 
B — B1/2, the Hall resistance g x is determined by the 
cyclotron frequency Wtot = Wc + ^1/2 f° r the: total mag- 
netic field B, with W1/2 = eBi/2/m* CFl due to the in- 
clusion of the Chern-Simons electric field. Introducing 
the mean free path Acf = T tr hk F /m* CFl we obtain 
q2 f = hk^ F /(e 2 n\ C F) and g\ y = B/(ne), so that the 
resistivity tensor is independent of tuq Fi and the Hall re- 
sistance is the same as that of the 2DES. Comparing with 
the corresponding Drude formulas for the 2DES at low 
B, we see that the ratio of the longitudinal resistances, 
Qq F J gf = y^Aei/AcF, directly reflects the ratio of the 
corresponding mean free paths. 

The macroscopic resistivity tensor g of the modulated 
system, relates the spatial average (j(x)) of the current 
density j(x) to the driving field, g(j(x)) — E' ). Ex- 
ploiting the solution for the homogeneous system and 
the continuity equation, which implies that j x (x) for the 
modulated system is independent of x, it can be shown 
that g differs from g h only in its xx-component. With 
some formal but exact manipulations (similar to those of 
ref. M) one finds 

(&* - qT)/Q™ = G/[l - G/(l + (wtotTfe) 2 )]. (7) 



In this equation 

2Ttr I r dp 



G 



2tt 



g(x,ip)x(x,<p) 



(8) 



where g(x,ip) = v F (dV / dx) / (2E F F ) + u) m v y is of first 
order in the modulation strength ecF) and x{ x t f) is the 
solution of the modified Boltzmann equation 



£>x - c\x] = g{x, <p) + uj 1/2 v y (jx/ji 1 - 1) 



(9) 



The basic differences between these results for the CF 
system and those for an electron system M derive from 
inclusion of the Chern-Simons electric field in Eq. (0): 



It is responsible for the occurence of w to t oc B in the 
denominator of Eq. (0), instead of uj c oc B — B1/2, and 
the last term in the right hand side of Eq. (0). This 
last term can be expressed in terms of G, j x /Jx ~ 1 = 
— G/(l + (wtot7"tr) 2 ), and requires a self-consistent solu- 
tion of Eqs. (Q)-(||). Inserting r tr = A C f/vf into Eq. (||) 
makes G and xa; independent of m CF , as it should be. 

Numerical solution of Eqs. (0)-(||) for the CF system 
and of the corresponding equations of Ref. M for the elec- 
tron system allows to calculate the resistance of the mod- 
ulated sample for strong applied fields B w -B1/2 and in 
the low B regime, respectively. To compare with experi- 
ment, we need to choose a reasonable set of parameters. 
Whereas in both B regimes the average density h and the 
modulation period a are the same, and also the modula- 
tion strengths 6cf ~ £ei should be similar, the parameters 
characterizing the collision operator (in our model r, b 
and p) are expected to be different. To reproduce the 
characteristics of the Weiss oscillations observed at small 
B J6[, we have to assume a very anisotropic differential 
cross section with a sharp peak in forward direction and 
a zero- 5 resistance of only a few Q (A c i ~ 50 /xm). This 
is reasonable since these high-mobility samples have a 
spacer thickness of about 60 nm, so that the donors pro- 
duce a smooth random potential in the plane of the 2DES 
that scatters the electrons predominantly under small an- 
gles. The CFs, on the other hand, are scattered not only 
by these small-amplitude donor-induced random poten- 
tial fluctuation, but also - and predominantly - by the 
large-amplitude random fluctuations of the effective mag- 
netic field that originates from the donor-induced density 
fluctuation. This results in Acf *C A i and a less pro- 
nounced forward scattering for the CFs. 

Figure [fl shows resistivity data for a sample with den- 
sity n = 1.82 ■ 10 11 cm -2 and period a = 400 nm near 
filling factor v = 1/2 w\- The thick short-dashed line is 
the experiment of Ref. |6| . The thin dash-dotted line is 
the theoretical curve using the approximations of Ref. [|S| , 
based on the Becnakkcr-type approximation (BA) men- 
tioned above. It neglects the modulation effects in the 
drift operator T> of the LBE, that causes a resistance cor- 
rection quadratic in ecf- Furthermore, the direct effect of 
the electric modulation on the CFs was ignored in this fit. 
Neglecting this direct effect too and using the BA, we can 
closely reproduce this symmetric curve of Ref. 8 using 
the same values for n, Qq, £cf, and using b = 0, p = 6 
for our model of the scattering cross section, Eqs. (|5|) 
and (0). If we include the direct electric modulation in 
addition to the magnetic one, we obtain the asymmetric 
thin dotted line. For small \B c fi \, this again closely repro- 
duces the results of Ref. M, but at large \B e s | the curve 
shows a distinct asymmetry, similar to the experimental 
curve, and similar to the asymmetry obtained for mixed 
electric and magnetic modulations in 2DESs at low B 
0,0. Solving Eqs. (R)-(||) for the same set of param- 
eters but without any approximation of the drift opera- 



tor, we obtain the thin solid curve of Fig. [j] This shows 
that approximating the drift operator T> by the one of 
the homogeneous system is insufficient for |-B e ff| < 0.5 T. 
Physically the BA means that the forces due to the mod- 
ulation fields are neglected against the Lorentz force due 
to the average magnetic field -B c ff- This omits the ef- 
fect of channeled orbits, i.e. of "snake orbits" |13[ along 
the lines of vanishing B c k{x) which, similar to the "open 
orbits" [J14J in the case of a pure electric modulation, 
cause a pronounced positive low-field magneto-resistance 
@0|. Since B 1/2 = 15 T ||, the effective magnetic mod- 
ulation has an amplitude ecF-Bi/2 ~ 0.45 T, so that the 
effect of channeled trajectories should be important for 
|-E?eff J < 0.45 T, and a second order approximation in ecf 
is not adequate for the B e g values shown in Fig. [j] 

Comparing the thin solid curve with the experimental 
one suggests that the pronounced V-shaped resistance 
minimum at i? ff = seen in the experiment is indeed due 
to CFs on channeled orbits. To achieve a better quan- 
titative agreement, we choose more realistic parameters. 
First we take b — 0.1 and p — 2, since then our model 
(pf) approximates closely (apart from the ip = diver- 
gence) the differential cross section derived by Aronov et 
al. |l5[ for random magnetic field scattering. This choice 
allows a good fit of the resistance near v = 1/2 for all 
samples available to us, indicating the same scattering 
mechanism in those samples. Then, assuming a density 
modulation of about 3%, we have to take Qq = 650 il 
(Acf ~ 1.3a), i.e. a larger value than the resistance of 
the unmodulated reference sample (qq = 270 fi) taken 
in Ref. M. This yields the thick solid curve from the 
exact solution of the LBE, which reproduces nicely all 
features of the experimental data. The larger g^ f value 
is reasonable since the etching procedure, in addition to 
the intended modulation, introduces unintended defects 
which increase the resistance. 

For the sample discussed so far (a — 400 nm), the 
mean free path of the CFs is only marginally larger than 
the modulation period, Acf ~ 1.3a, so that commensu- 
rability effects are visible only as weak shoulders (thick 
lines of Fig. [j]). To see these effects more clearly, we have 
investigated a sample with a smaller modulation period, 
a = 275 nm, so that we expect Acf ~ 1.9a. Experimen- 
tal |lj] and theoretical results (with the same anisotropy 
parameters b and p as in Fig. ||) are plotted in Fig. g. 
Indeed, the commensurability effects near |_B G ff| = 0.5 T 
are more pronounced. Again, the low field resistance (in- 
set of Fig. 0) and that near v = 1/2 can be reproduced 
theoretically with a consistent set of parameters. 

In conclusion, our numerical solution of the linearized 
Boltzmann equation for composite Fermions without the 
Beenakker-type approximation used in previous work 
P,|l, and thus the inclusion of the effect of "channeled 
orbits" , reproduces the typical V-shaped resistance min- 
imum at v — 1/2 (S c ff = 0) and yields good overall 
agreement with experiment, consistent with the low-field 



magnetoresistance oscillations. Including further the su- 
perposition of the original electric superlattice and the 
induced effective magnetic superlattice, we can also re- 
produce the asymmetry of the resistance curve around 
B cS = 0. 

We are grateful to J. Smet for providing the unpub- 
lished data of Fig. and for helpful discussions. This 
work was supported by BMBF Grant No. 01BM622. 



FIG. 2. Resistivity near filling factor 1/2, a — 275 nm, 
n = 1.98 • 10 11 cm~ 2 ; thick dashed line: experimental data, 
thick solid line: calculation for Qq = 700 fi, and ecf = 3.8%, 
cross section parameters p — 2 and b = 0.1, as in Fig. m 
(thick lines), thin dotted line: Beenakker-type approximation 
(same parameters) . The inset shows the low-field magnetore- 
sistivity; dashed line: experiment, solid line: calculation for 
go = 7.5 n, eel = 3.6%, b = 0, and p = 15. 
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FIG. 1. Resistivity near filling factor 1/2, a = 400 nm, 
n — 1.82 ■ 10 11 cm -2 ; thick dashed line: experiment of Ref. 
[pj, thin dash-dotted line: theory of Ref. |8j (parameters: 
Jo = 270 Q and ecf = 2.6%), thin solid line: results of 
full calculation (including channeled orbits) with p — 6 and 
6 = 0, thin dotted line: Beenakker approximation for same 
parameters, thick solid line: full calculation for Qq f = 650 Q., 
e CP = 3.5%, p = 2 and b = 0.1. The inset shows the low-field 
magnetoresistance; dashed line: experiment, solid line: full 
calculation based on Ref. M for Qq = 4.5 ft and e e i = 2.8%, 
6 = and p — 15. 
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